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This  report  describes  a  study  to  identify  characteristics  which 
induce  thermal  stratification  in  liquid  thermal  storage,  and  to  evaluate 
solar  energy  system  performance  as  a  function  of  the  degree  of  stratifi¬ 
cation. 

It  was  determined  that  for  efficient  use  of  thermal  stratification 
it  was  necessary  to  (1)  introduce  hot  fluid  at  the  top  of  the  liqu 
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^  storage  tank  and  to  add  cold  fluid  at  the  bottom/of  the  tank,  (2)  with  a 
cylindrical  tank,  have  a  length/diameter  ratio($)2.0,  (3)  use  mathemati¬ 
cal  correlations  to  determine  allowable  fluid  inlet  velocities  and  tem¬ 
peratures,  and  (4)  use  storage  tank  material  that  has  a  thermal  conduc¬ 
tivity  less  than  that  of  the  storage  fluid. 

The  mathematical  correlations  described  in  this  report  allow  stra¬ 
tification  occurrence  to  be  predicted  and  can  be  used  to  estimate  the 
sharpness  of  the  thermocline  based  on  tank  inlet  and  outlet  conditions, 
fluid  properties,  and  storage  tank  geometry. 
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INVESTIGATION  OF  METHODS  TO  PREDICT  THERMAL 
STRATIFICATION  AND  ITS  EFFECT  ON  SOLAR  ENERGY 
SYSTEM  PERFORMANCE 


1  INTRODUCTION 


Background 

The  high  costs  of  conventional  and  alternate  energy  sources  require 
that  they  be  used  efficiently;  one  way  of  insuring  energy  efficiency  and 
conservation  is  to  store  excess  energy  for  later  use.  Energy  storage  is 
sometimes  a  matter  of  physical  necessity,  as  in  solar  energy  systems, 
and  sometimes  a  matter  of  economic  necessity,  as  when  electric  utilities 
offer  a  peak  billing  program. 

Because  much  of  the  energy  consumed  today  exists  at  one  time  or 
another  as  thermal  energy,  there  is  an  increasing  use  of  and  need  for 
thermal  energy  storage. 

Thermal  energy  may  be  stored  in  various  ways:  (1)  as  sensible  heat 
by  increasing  of  a  material's  temperature,  (2)  as  latent  heat  by  chang¬ 
ing  the  phase  of  a  material  (i.e.,  melting,  boiling),  or  (3)  as  a  combi¬ 
nation  of  (1)  and  (2). 

Because  sensible  heat  storage  is  simple  and  inexpensive,  it  is 
currently  the  most  common  energy  storage  method;  the  materials  most 
often  used  to  store  energy  as  sensible  heat  are  water  and  rock.  Sensi¬ 
ble  storage  also  sometimes  allows  temperature  stratification,  which  may 
improve  the  performance  of  systems  using  the  sensible  heat  storage 
met hod. * 


Objective 


The  objective  of 
induce  stratification 
system  performance  as 


this  report  is  to  identify  characteristics  which 
in  liquid  thermal  storage,  and  to  evaluate  energy 
a  function  of  the  degree  of  stratification. 


1  T.  D.  Brumleve,  Sensible  Heat  Storage  ^  j-i quids.  Plowshare  and  Trans¬ 
ducer  Technology  Division  8184,  Report  SLL-/3-0263  (Sandia  Labora¬ 
tories,  July  1974). 
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Approach 

An  experimental  apparatus  was  designed  and  constructed  to  evaluate 
the  impact  of  inlet/outlet  conditions,  fluid  properties,  and  tank 
geometry  on  the  fluid  temperature  profile.  Experiments  and  analyses 
were  performed  based  on  correlating  experimental  data  with  previously 
identified  dimensionless  parameters.2  In  addition,  several  literature 
searches  were  undertaken.  Based  on  the  results  of  the  experiments  and 
the  available  literature,  a  computer  simulation  of  a  solar  energy  system 
was  performed. 


2  A.  Cabelli,  "Storage  Tanks,  a  Numerical  Experiment,"  Solar  Energy,  Vol 
19,  No.  1  (Pergamon  Press,  Inc.,  1977),  pp  45-54. 
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2  CHARACTERISTICS  WHICH  INDUCE  STRATIFICATION 


Definition  of  Thermal  Stratification 

Thermal  stratification  exists  in  storage  fluid  when  there  is  a  tem¬ 
perature  differential  between  top  and  bottom  fluids.  This  temperature 
differential  is  caused  by  buoyant  forces  resulting  from  the  density 
difference  between  hot  and  cold  fluid. 

This  definition  is  quite  broad  and  would  include  all  the  conditions 
shown  in  Figure  1.  In  fact,  the  buoyant  forces  are  so  pervasive  that 
some  degree  of  stratification  exists  in  almost  all  storage,  except  for 
storage  with  very  large  mixing  forces. 

Figure  la  shows  the  components  of  stratified  storage  as  it  was 
defined  for  this  study.  There  is  a  hot  mixed  layer,  a  cold  mixed  layer, 
and  a  boundary  region  (thermocline)  separating  the  hot  and  cold  mixed 
layers.  The  lower  limit  of  stratification  is  reached  when  the  thermo¬ 
cline  occupies  the  entire  tank  (Figure  lb).  Figures  lc  and  Id  show  par¬ 
tially  and  fully  mixed  conditions. 


Characteristics  Which  Induce  Stratification 


The  phenomenon  of  thermal  stratification  in  thermal  storage  is  com¬ 
plex.  The  important  elements  of  the  thermal  stratification  problem  are: 
buoyant  forces,  mixing  forces,  thermal  diffusivity,  thermal  conduc¬ 
tivity,  inlet  configurations,  and  tank  geometries.  Also  important  are 
the  variations  of  some  of  these  elements  with  time  as  storage  is  used  in 
a  system. 

The  results  of  the  research  to  date  indicate  that,  in  general,  the 
following  elements  enhance  stratification: 

1.  Warm  fluid  inlet  at  top,  cold  fluid  inlet  at  bottom 

2.  Increasing  the  temperature  differential  between  tank  fluid  and 
incoming  fluid 

3.  Reducing  the  inlet  fluid  velocity 

4.  Decreasing  the  tank  wall  conductivity 

5.  Increasing  the  L/D  ratio 

where: 

L  =  length  of  cylindrical  tank  (m) 

D  =  diameter  of  cylindrical  tank  (m). 
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Figure  1.  Examples  of  temperature  profiles. 


Effects  of  Inlet  Temperatures,  Velocities,  and  Tank  L/D  Ratio 

There  are  twc  correlations  available  which  predict  the  occurrence 
of  storage  stratification  as  a  function  of  operating  temperatures,  inlet 
velocities,  and  tank  dimensions.  The  correlations  were  developed  con¬ 
currently  by  different  methods,  but  they  are  very  similar.  Both  of  the 
correlations  are  based  on  experimental  data  gathered  under  the  following 
conditions: 

1.  Constant  AT  between  inlet  fluid  and  tank  fluid 

2.  Constant  flow  rate 

3.  Well-insulated  tanks  with  low  tank-wall  conductivity. 

The  first  of  these  correlations  is  that  proposed  by  Lavan  and 
Thompson. ^  in  their  experiments,  a  cylindrical  plexiglass  tank  was 
filled  first  with  hot  water  at  105°F  (40.6°C),  then  with  cold  water  at 
60°F  (15°C),  which  was  pumped  into  the  tank  through  a  horizontal  inlet 
at  the  bottom  of  the  tank.  This  corresponds  to  the  condition  of 
discharging  energy  from  the  tank.  As  cold  fluid  displaced  hot  fluid  in 
the  tank,  the  exit  fluid  temperature  was  measured.  To  quantify  their 
results,  Lavan  and  Thompson  introduced  an  extraction  efficiency  defined 
as: 


n  =  Q-t/V 


[Eq  1] 


(where: 

n  is  the  extraction  efficiency  (dimensionless) 

Q  is  the  volumetric  flow  rate 
V  is  the  tank  volume 

t  is  the  time  it  takes  the  initial  inlet/exit  temperature 
difference  as  measured  at  the  exit  to  drop  some  preassigned 
value. 

In  these  experiments,  a  10  percent  drop  was  assumed,  i.e.: 

At  time  t,  (Z^ — Jj.n  1*1*1)  =  o.9  [Eq  2] 

1  in  '  initial 

I  Various  experiments  were  performed  with  different  values  of  L/D,  AT,  and 

l  inlet  velocity.  The  results  are  summarized  in  Figure  2. 

1 

\ 


3  Z.  Lavan  and  J.  Thompson,  "Experimental  Study  of  Thermally  Stratified 
Hot  Water  Tanks,"  Solar  Energy,  Vol  19,  No.  5  (1977),  pp  519-524. 
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The  data  points  indicated  by  A  in  Figure  2  were  taken  using  a  hor¬ 
izontal  distributor  constructed  of  0.7  in.  (17.7  mm)  diameter  copper 
tubing  with  fifty  0.112-in.-  (2.84-mm)-diameter  holes  drilled  in  two 
lines  along  the  tube  on  a  0.25-in.  (6.35-mm)  center.  Figure  2  shows  the 
improvement  in  extraction  efficiency  obtained  with  the  distributor.  A 
second  correlation  was  developed  based  on  data  from  a  previous  study4 
and  data  gathered  during  this  study.  This  correlation  was  developed  by 
observing  the  temperature  profile  along  the  center  line  of  a  cylindrical 
fiberglass  tank  with  an  L/D  of  2.  The  inlets  and  outlets  were  horizon¬ 
tal  and  made  of  1/4-  and  1/2-in.  (6.35-  and  12.7-mm)  schedule  40  steel 
pipe.  Experiments  were  performed  under  conditions  of  charging  and 
discharging.  In  the  charging  experiments,  the  tank  was  filled  initially 
with  water  ranging  in  temperature  from  70  to  150°F  (21  to  65°C).  Warmer 
fluid  was  then  pumped  in  at  the  top  of  the  tank  as  cooler  fluid  was 
removed  at  the  bottom.  The  differential  between  inlet  fluid  temperature 
and  initial  tank  temperature  ranged  from  10  to  50°F  (5.56  to  27.78°C). 
Flow  rates  ranged  from  0.2  to  1.5  gpm  (0.76  to  5.68  £/min). 

In  the  discharging  experiments,  the  initial  tank  temperature  ranged 
from  80  to  135°F  (26  to  57°C).  In  the  discharging  experiments,  cooler 
fluid  was  pumped  in  at  the  bottom  of  the  tank.  The  differential  between 
the  inlet  and  initial  tank  fluid  temperature  ranged  from  10  to  70°F 
(5.56  to  38.89°C).  The  apparatus  used  in  these  experiments  is  shown  in 
Figure  3. 


Figure  3.  Experimental  apparatus. 


4  B.  Sliwinski,  A.  Mech,  and  T.  S.  Shih,  "Stratification  in  Thermal 
Storage  During  Charging,"  Proceedings,  Sixth  International  Heat 
Transfer  Conference,  Vol  4  (Toronto,  Canada,  1978),  pp  149-154. 


The  temperature  profiles  observed  in  these  tanks  indicated  that  the 
depth  at  which  the  thermocline  formed  varied  from  experiment  to  experi¬ 
ment.  Therefore,  the  depths  at  which  thermocl ines  formed  were  used  to 
quantify  the  occurrence  of  stratification.  Figure  4  shows  some  of  the 
profile  types  observed  during  the  experiments.  The  dimensionless  mixing 
pa rameter 


L*/L  [Eq  3] 

is  introduced  where  L*  is  the  depth  where  the  thermocline  forms  and  L  is 
the  distance  between  the  inlet  and  outlet. 


The  parameter  L*/L  was  correlated  with  a  dimensionless  parameter 
known  as  the  Richardson  number.  This  number  is  a  ratio  of  buoyant  to 
mixing  forces: 


Ri 


g  pat0l 


[Eq  4] 


where: 

Ri  =  Richardson  number  _ 

g  =  acceleration  of  gravity  (m/s^) 

e  =  volume  expansivity  (l/°C) 

v.j  =  inlet  velocity  (m/s) 

aTq  =  difference  between  inlet  fluid  and  initial 

tank  fluid  temperature  (°C). 

This  correlation  is  shown  in  Figure  5. 

When  L*/L  =  1.0,  the  tank  is  fully  mixed;  when  L*/L  =  0.0,  stratif¬ 
ication  occurs  immediately  at  the  inlet  of  the  tank. 


Limits  to  Present  Correlations 


The  correlations  described  previously  improve  the  designer's  abil¬ 
ity  to  predict  stratification,  but  they  have  limitations.  The  most 
important  limitation  is  the  inability  to  predict  the  effect  of  varying 
input  temperature  on  the  fluid  temperature  profile.  Preliminary  experi¬ 
ments  were  performed  to  study  the  effect  of  varying  inlet  temperatures. 
The  results  of  these  preliminary  experiments  are  shown  in  Figure  6. 
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Figure  4.  Examples  of  the  types  of  profiles  observed  experimentally. 


Discharge  Data 
Charge  Data 

t#  -»/* 

=  0.0313  Ri  for  0001  <  Ri  <  10.0 


Figure  5.  Occurrence  of  stratification  as  function  of  Ri . 


Sharpness  of  Thermocline 

The  sharpness  of  the  thermocline  characterizes  the  degree  of  stra¬ 
tification.  The  sharpness  can  be  used  to  predict  thermal  tank  material 
stresses.  The  sharpness  of  thermocline  is  given  by 


aT/aZ/ ATq/L 


[Eq  5] 


where: 


aZ  is  the  thickness  of  the  linear  portion 
of  the  thermocline 

AT  is  the  temperature  difference  across  the 
linear  portion  of  the  thermocline  (Figure  7). 

Eq  5,  correlated  with  the  Richardson  number  is: 

Ri  -  96AT0L 

Rid  "  — jr- 

u 


[Eq  6 


LJ 


and  with  the  Peclet  number  defined  as: 


Pe  «  Re^Pr 


[Eq  7] 


where : 

ReD  is  the  Reynolds  number  based  on  the  tank  diameter 
Pruis  the  Prandtl  number 

The  Peclet  number  is  the  ratio  of  heat  transfer  by  convection  to 
heat  transfer  by  conduction.  These  correlations  are  shown  in  Figures  8 
and  9. 


At  this  time  it  is  difficult,  based  on  the  scantiness  of  the  data, 
to  place  a  high  degree  of  confidence  in  these  correlations.  However,  as 
a  guide  to  estimating  the  sharpness  of  the  thermocline,  aT/aZ/aT /L 
ranged  from  8  to  3,  with  an  average  from  all  experiments  of  4.09,  and  a 
standard  deviation  of  1.23. 


Effect  of  Tank  Wall  Conductivity  on 
Thermal  Stratification 


The  decay  of  the  thermocline  is  a  function  of  the  conduction 
between  the  hot  and  cold  layers  through  the  fluid  and  through  the  tank 
wall.  When  the  tank  wall  has  a  thermal  conductivity  higher  than  that  of 
the  fluid,  the  wall  conduction  effects  can  be  6  to  7  times  larger  than 
fluid  conduction.5  In  general,  tank  wall  material  should  be  selected  to 
have  a  thermal  conductivity  less  than  that  of  the  storage  fluid. 


5 


C.  W.  Miller,  The  Effect  of  a  Conducting  Wall  on  a  Stratified  Fluid 
a  Cylinder,  Aerospace  Industries  Association  of  America  (AlAA), 
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3  THE  EFFECT  OF  THERMAL  STRATIFICATION  ON  THERMAL 
SYSTEM  PERFORMANCE 


General 


The  effect  of  stratification  on  thermal  system  performance  is  the 
subject  of  current  research.  The  literature  search  conducted  during 
this  study  indicated  stratification  improves  system  performance,  but 
discovered  no  consensus  as  to  the  amount  of  improvement. 


Mechanisms 


Thermal  storage  stratification  affects  system  performance  through 
two  mechanisms: 

1.  Higher  charging  and  discharging  rates 

2.  Constant  outlet  temperatures. 

The  following  equations  show  how  stratified  storage  allows  for 
higher  charging  and  discharging  rates  than  mixed  storage;  i.e.,  for  a 
mixed  tank  with  constant  temperature  input  and  constant  flow  rate,  the 
energy  content  as  a  function  of  time  is 


-mt 

E  =  aTq  MCp  (1  -  e  M  ) 


[Eq  8] 


where: 

E  =  energy  content  as  a  function  of  time 
t  =  time 

Cp  =  specific  heat 
m  =  mass  flow  rate  (kg/s) 

M  =  mass  of  storage  fluid  (kg) 


For  a  perfectly  stratified  tank,  the  energy  content  as  a  function  of 
time  is 


E  =  AT0mCpt 
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[Eq  91 


The  ratio  E/E  (where  E  is  the  minimum  energy  capacity  of  the 
storage)  is  shown  as  a  functiorrof  the  tank  time  constant  M/m  in  Figure 
10.  Figure  10a  shows  that  the  stratified  tank  is  fully  charged  in  one 
time  constant,  whereas  it  takes  about  4  to  5  time  constants  to  fully 
charge  the  mixed  tank. 


The  equation  for  discharging  a  mixed  tank  is: 


-mt 


E  =  AT0MCpe  " 

The  equation  for  discharging  a  stratified  tank  is: 


E  =  aTq  Cp  (M  -  mt) 


[Eq  10] 


[Eq  11] 


Eqs  9  and  11  are  plotted  the  same  way  in  Figure  10b,  which  shows  that 
the  stratified  tank  is  fully  discharged  in  one  time  constant,  whereas  it 
takes  about  4  to  5  time  constants  to  fully  discharge  the  mixed  tank. 


The  rapid  charge  and  discharge  characteristics  of  stratified  ther¬ 
mal  storage  will  dynamically  affect  a  system,  but  the  degree  of  effect 
will  vary  from  system  to  system.  In  general,  there  will  be  a  decrease 
in  the  operating  time  of  pumps  and  fans  used  to  add  or  remove  energy 
from  storage.  Also,  the  storage  will  be  better  suited  to  supply  systems 
requiring  high  energy  input  rates. 


The  constant  outlet  temperature  of  stratified  storage  also  affects 
system  performance.  Figure  11  shows  the  outlet  temperatures  of  mixed 
and  stratified  tanks  during  charge  and  discharge.  Note  that  the  outlet 
temperature  of  the  stratified  tank  remains  constant  through  nearly  the 
entire  charge  or  discharge  cycle,  but  the  mixed  tank  outlet  temperature 
varies. 


Constant  outlet  temperature  affects  different  systems  in  various 
ways,  but,  in  general,  when  a  system  requires  low  supply  temperatures 
(e.g.,  a  solar  collector)  or  constant  high  temperatures  (e.g.,  an 
absorption  chiller),  stratified  storage  is  better  than  mixed  storage. 


Effects  on  Solar  Energy  System  Performance 

The  .iterature  search  conducted  during  this  study  indicated  that 
research  into  system  performance  usually  was  done  either  by  experiment¬ 
ing  with  actual  solar  energy  systems  or  by  using  computer  simulations. 
Both  approaches  have  encountered  problems. 

Many  actual  system  experiments,  because  of  lack  of  understanding  of 
the  stratification  problem,  did  not  properly  use  the  stratified  storage 
energy  system.  In  the  case  of  computer  simulations,  the  one-dimensional 
models  generally  used  were  not  precise  enough  to  model  thermal  storage. 
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Charging 


Discharging 


Figure  10.  Charging  and  discharging  rates  of  mixed  and  stratified  storage 


Charge  Time 


Figure  11.  Output  temperatures  of  mixed  and  stratified 
storage  during  charge  and  discharge. 


Computer  simulations  using  two-  and  three-dimensional  models,  while 
becoming  very  accurate,  are  not  suited  for  detailed  system  simulation.6 

Some  results  of  recent  research  are  summarized  in  Table  1.  Table  1 
shows  that  the  amount  of  improvement  in  system  performance  varied  from  0 
to  more  than  30  percent.  Figure  12  shows  the  improvement  in  system  per¬ 
formance  calculated  using  a  modified  TRNSYS?  model  with  a  lengthened 
time  step. 

In  this  simulation,  the  annual  performance  of  a  solar  heating  and 
cooling  system  was  projected. 

The  storage  tank  model  used  was  a  six-segment  tank  similar  to  that 
used  in  past  TRNSYS  simulations  with  the  exception  of  the  use  of  a 
lengthened  time  step  to  calculate  heat  flow  between  each  segment. 

The  lengthened  time  step  was  used  based  on  comparisons  of  tank 
fluid  temperature  profiles  predicted  by  the  model  for  charging  and 
discharging  and  experimental  data  from  the  same  conditions.  The  best 
agreement  was  found  when  a  lengthened  time  step  was  used. 


6 

7 


Personal  communication  between  B.  J.  Sliwinski  of  the  U.S.  Army  Con¬ 
struction  Engineering  Research  Laboratory  (CERL)  and  E.  I.  H.  Lin  of 
the  Argonne  National  Laboratory,  August  7,  1979. 

"TRNSYS  --  A  Transient  Simulation  Program,"  American  Society  of  Heat¬ 
ing,  Refrigeration,  and  Air-Conditioning  Engineers  (ASHRAE)  Transac¬ 


tions,  Vo!  82.  Part  1  (19761. 
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Figure  12.  Improvement  in  performance  of  sii 
solar  heating  and  cooling  system 


Table  1 


Results  of  Research  to  Determine  Effect  of 
Stratification  on  Solar  Energy  System  Performance 


Type  of 


Researcher 

Study 

Result 

Reference 

Namkoong 

Experimental 

No  improvement 

* 

S.  T.  Wu 

1-D  TRNSYS 
simulation 

6%  improvement  in 
collector  efficiency 

** 

Y.  NakaJima 
(et  al . ) 

Simulation  and 

experimental 

system 

20-30%  improvement 
in  collector  efficiency 

★  ** 

Van  Koppen 

Experimental 

system 

15%  system  improvement 

+ 

E.  H.  Lin 

3-D  simulation 
analysis 

30%  system  improvement 

+  + 

Walton  & 

SI iwinski 

1-D  TRNSYS 
simulation 
with  lengthened 
time  step 

15-20%  system 
improvement 

Results  of  the  study 
documented  in  this 
report 

*  D.  Namkoong,  Temperature  Distribution  of  a  Hot  Hater  Storage  Tank 
in  a  Simulated  Solar  Heating  and  Cooling  System,  NASA  Report 
TM  X-73549  (National  Aeronautics  and  Space  Administration  [NASA], 

Lewis  Research  Center,  November  1976),  p  7. 

**  S.  T.  Wu,  The  Theoretical  and  Experimental  Study  of  Liquid  Storage 
Tank  Thermal  Stratification  for  Solar  Energy  Systems,  DOE  Report  No. 
COO/4479-1  (Department  of  Energy  LDOEJ,  February  1978),  p  1. 

***  Y.  NakaJima,  "Design  and  Performance  of  Thermal  Storage  water  Tanks," 
Heliotechnique  and  Development,  Proceedings  of  International 
Conference,  Vol  1,  Dhahran,  Saudi  Arabia  (November  1975),  pp  508-510. 

+  C.  W.  J.  Van  Koppen,  L.  S.  Fischer,  and  A.  Dijkmans,  "Stratification 
Effects  in  the  Long  and  Short  Term  Storage  of  Solar  Heat,"  Proceedings, 
1978  meeting  of  Amercian  Sec.  International  Solar  Energy  Society 
(ISES)  (Denver,  CO,  1978). 

++  E.  T.  H.  Lin,  On  Thermal  Energy  Storage  Efficiency  and  the  Use  of 
COMMIX-SA  for  Tts  Evaluation  and  Enhancement,  Solar  Energy  Storage 
Options  Workshop,  San  Antonio,  TX,  1979,  pp  7-8. 
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4  CONCLUSIONS  AND  RECOMMENDATIONS 

To  assure  efficient  use  of  thermal  stratification  in  liquid  thermal 
storage,  it  is  necessary  to: 

1.  Introduce  hot  fluid  at  the  top  of  the  liquid  storage  tank  and 
add  cold  fluid  at  the  bottom  of  the  tank. 

2.  With  a  cylindrical  tank,  use  an  L/D  ratio  2  2.0. 

3.  Determine  allowable  fluid  inlet  velocities  and  temperatures 
from  correlations  in  Figure  2  or  Figure  5. 

4.  Use  tank  material  that  has  a  thermal  conductivity  less  than 
that  of  the  storage  fluid. 

Information  about  thermal  stratification  in  liquid  sensible  storage 
is  incomplete,  but  rapidly  growing.  Correlations  have  been  developed 
which  allow  the  designer  to  predict  stratification  occurrence  and  esti¬ 
mate  the  sharpness  of  the  thermocline  based  on  inlet  and  outlet  condi¬ 
tions,  fluid  properties,  and  tank  geometry.  These  correlations  are  not 
limited  in  application  to  solar  energy  systems.  On  the  contrary,  they 
are  most  applicable  to  conventional  systems  which  have  operating  condi¬ 
tions  similar  to  the  experimental  conditions  from  which  the  correlations 
were  derived,  i.e.,  constant  input  temperature  and  flow  rate.  They  can 
also  be  applied  to  solar  energy  systems  with  the  understanding  that  in 
solar  energy  systems,  the  tank  input  temperature  from  the  collectors  may 
not  be  constant. 

The  results  of  CERL's  simulation  and  work  in  the  literature  indi¬ 
cate  that  the  effects  of  strati fication  on  the  performance  of  solar 
energy  systems  are  positive,  showing  an  improvement  in  system  perfor¬ 
mance  of  between  20  and  30  percent  compared  to  a  mixed  tank  system. 

It  is  recommended  that: 

1.  Systems  using  sensible  heat  storage  should  be  designed  to  max¬ 
imize  stratification  occurrence. 

2.  The  response  of  stratified  storage  to  time  varying  inlet  tem¬ 
peratures  which  fluctuate  above  and  below  the  temperature  of  the  mixed 
region  in  the  tank  should  be  determined.  In  addition,  there  is  urgent 
need  for  an  accurate,  quasi  one-dimensional  computer  simulation  model  of 
stratification  that  can  be  used  in  system  simulations  to  assess  the 
impact  of  using  stratification  on  the  first  cost  of  other  system  com¬ 
ponents.  It  is  also  important  to  determine  suitable  tank  construction 
materials  to  insure  low  tank  wall  thermal  conductivity. 


..I.-  - - -V:.  ■  ■  -»*■  -i.  1  H 
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NOMENCLATURE 


L  =  Length  of  cylindrical  tank  (m) 

D  =  Diameter  of  cylindrical  tank  (m) 

T  =  Temperature  (°C) 

n  =  Extraction  efficiency  (dimensionless) 

o 

g  =  Acceleration  of  gravity  (m/s  ) 

8  =  Volume  expansivity  (1/°C) 

Vi  =  Inlet  velocity  (m/s) 

aT  =  Difference  between  inlet  fluid  and  (°C) 
initial  tank  fluid  temperature 

aT  =  Temperature  difference  across  thermocline  (°C) 

AZ  =  Thickness  of  the  linear  portion  of  the  thermocline 

Cp  =  Specific  heat  (i2Hlg£) 

•  ^ 
m  =  Mass  flow  rate  (kg/s) 

M  =  Mass  of  storage  fluid  (kg) 

Rej)  =  Reynolds  number  based  on  tank  diameter 
bulk  fluid  velocity  and  average  of 
hot  and  cold  fluid  temperature  (dimensionless) 

Pr  =  Prandtl  number  based  on  average  of  hot  and 
cold  fluid  temperature  (dimensionless) 

Vp  =  Bulk  fluid  velocity  based  on  tank  diameter  (m/s) 
Red  =  Inlet  Reynolds  number 
Grp  -  Tank  Grashof  number 
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